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Abstract 

In this paper, a new variational formulation based on discontinuous Galerkin 
technique for a reaction-diffusion problem is introduced, and the discontinu- 
ous Galerkin technique of this work is different from the general discontinuous 
Galerkin methods. The well posedness of the new formulation is given. Fi- 
nally, it is pointed that the new variational formulation will be helpful to 
design better hybrid numerical methods which will not only strongly stable 
in spatial variable and absolutely stable in temporal variable but also be 
optimally convergent. 

Keywords: Discontinuous Galerkin technique, variational formulation, 
inf-sup condition. 

1. Introduction 

In this work, we propose a new variational formulation based on discon- 
tinuous Galerkin technique for a reaction-diffusion problem within a new 
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function space setting, the reaction-diffusion problem is 



V ■ (K(x)Vu) +u 



f, in fi, 



(1.1) 



u 



0, on dfl, 



(1.2) 



where / is a real- valued function in L 2 (Q) and < K < K(x) < K\. 

The problem fll.ip -f TOj) is an important and basic mathematical model, 



was introduced by [2j, and since that time there has been an active de- 
velopment of DG methods for hyperbolic and nearly hyperbolic problems, 
resulting in a variety of different methods. For elliptic and parabolic equa- 
tions, discontinuous finite elements were proposed by many researchers, such 



The key idea of the paper is to propose a new variational formulation 
based on discontinuous Galerkin technique, and the discontinuous Galerkin 
technique is different from the general discontinuous Galerkin methods. The 
formulation satisfies a local conservation property, and we prove well posed- 
ness of the new formulation by proving and using inf-sup condition. 

The paper is organized as follows. In Section 2, we introduce the new 
weak formulation of the problem f ll.ip - fll.2p . In Section 3, we investigate 
the well posedness of the variational formulation, which includes the conti- 
nuity property of the bilinear form and the inf-sup condition. Finally, some 
concluding remarks are summarized. 




as [3r|l6|. 
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2. A new variational formulation 

Let Q C M 2 be a bounded open domain with Lipschitz boundary dQ and 
let {Ph} be a family of regular partitions of Q into open elements E such 

that Q = int( [j E). 

EeP h 

The following notations will be used in our further considerations. Denote 
h = max^gp^ where kg = dia,m(E). The set of all edges of the partition 
P h is given by e h = {^ k }, k = 1, . . . , N edge , where N edge reprents the number 
of edges in the partition Ph. The interior interface Tint is then defined as the 
union of all common edges shared by elements of partition Ph, that is, 

-^edge 

Tint = |J le\OQ. 
k=l 

For the sake of clarity in the notation, the jump and average operators 
are defined by 

H = V^cdEi ~ v\ 7eCdEj , (v) = -{v^cdEi + v\ JeC dEj), % > j, (2.1) 

where 7 e = int(dEi D dEj) is the common edge in 2D (or interface in 3D) 
between two neighbouring elements, see Figure HJ 



3 




Figure 1: Geometrical definitions of neighbouring elements. 



First, we introduce the following broken Sobolev space: 



M{P h ) = {ve L 2 (Q)\v \ E e H(A, E),VE e P h , [Vv • n] e L 2 (r int )}, 



where 



H(A,E) = {ve L 2 (E)\V-Vv e L 2 (E)} C H\E). 



Notice here, that v e H(A, E) implies Vv ■ n e H 1 l 2 (dQ). The norm 
on M(Ph) is defined as 



-lll 2 = E{lkll^ + ^ll^)v.-iu||^ 1/2(aB) } 



E£P h 



+ a-\\lK(x)Vv»}\\l 2iri] 



p 



(2.2) 



where we denote that \\v\\ 2 = J E \K(x)\\Vv \ 2 dx+J E \v \ 2 dx = \\K2\7 



ll^lli 2 ^)' an< ^ one can eas ily prove that norms || ■ ||* and || ■ Hh 1 ^) are equiv- 
alent. The parameter p6R that is introduced here represents the minimum 
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of all of the local orders of polynomial approximations pe in the partition 
Ph- The parameters u, A, 6, ( are greater than or equal to zero and that the 
subsequent norms in (12. 2p are defined as 

\\ u \\H-v\dE) = su p — n — n ' ( 2 - 3 ) 

v eH 1 / 2 (dE) \\ L P\\m/ 2 (dE) 
IMItfVaps) = inf HU (2.4) 

J0"W=<P 

where (•, •)_ 1 / 2 xi/2 se denotes the duality pairing in H~ 1 ^ 2 (dE) x H^^dE), 
namely, 

(«,u)-i/2xi/2,as= / wfrfs. (2.5) 
Joe 

And 70 denotes the trace operator 

7o : H\E) ->• H l ' 2 {dE). 

Now, the choice for the space of test functions, V, is the completion of Ai(Ph) 
with respect to the norm 1 1 1 • 1 1 1 . 

The new discontinuous variational formulation, within this new function 
space setting, is then stated as follows: 

Find u E V, s.t, B(u, v) = L(v), \/v G V, (2.6) 
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where the bilinear form B(u,v) and linear form L(y) are defined as 
B(u,v) = ^2 { J (K(x)Vu-Vv + uv)dx 

(v(K(x) 



E€P h 



+ / ((v)[K(x)Vu ■ n] - (u)[K(x)Vv ■ n])ds 
Jr int 

r h x 

+ / a—[K(xWu-n][K(xWv -n}ds, (2.7) 
Jr int P c 

L{v) = [ fvdx. (2.8) 



Jn 

The formulation (JHHJ) (or VBVP (l2lj) ) is closelvrelated to the DG 
method formulation by Oden, Babuska and Baumann 17|. In fact, choosing 
the subspace V(Ph) of V of function with fluxes Vv ■ n e L 2 (dE), and using 
the following identities: 

^2 / v(K(x)Vu ■ v)ds 

E£P h ^ 9E 

= / [v(K(x)Vu ■ n)]ds + / v{K(x)Vu ■ n)ds, (2.9) 
Jr int ' Jan 

and 

[v{K{x)Vu ■ n)] = (K(x)Vu ■ n)[v] + (v)[K(x)Vu ■ n], (2.10) 



we can get the DG formulation of [17|. The only difference would then be 
;he addition of the last term in (12.71) . This term has been incorporated in 



18L Il9| . where it is accompanied by the jumps of the function [v] across the 
element interfaces. We replace the [v] jumps by the [Vv ■ /j] jumps, in order 
to prove both continuity and Inf-Sup properties of the bilinear form with 
respect to the space V, in which the norm is defined as 1 1 1 ■ 1 1 1 . 
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3. Well posedness of the new variational formulation 

In this section, we establish the well posedness of the variational formula- 
tion ( 12. 6p . Thus, we show that the solution of the problem (jl.ip - (ll.2p is also 
a solution to the weak problem. And we prove the existence, uniqueness and 
the continuous dependence on the input data of the solution to the varia- 
tional formulation ( 12. 6 p . Essential in some of these proofs are the continuity 
inf-sup conditions of the bilinear form (12. 7p . 



Now, we introduce an important lemma [16] as follows. 

Lemma 3.1. If u G H(A,Q), then u and (V« ■ n) are weakly continuous 
across the element interface T mt in the sense that 

I [u]ipds = Q, / [Vtt-n]^s = 0, (3.1) 

where ip belonges to D{^ e ), and 7 e = int(dEi G dEf) C r^t. 

Theorem 3.1. Let u be the solution of the problem / li.l)) -/ TOj) . Then u is a 
solution to the variational formulation 112.6]) as well. 

Proof. If we restrict (II. ip to an element E G Ph, multiply this local equation 
by a test function (pg G H 2 (E), integrate over the element E, and apply 
Green's identity, we get 

/ (K(x)Vu ■ VtpE + u<pE)dx — \ (Pe{K{x)Vu ■ fi)ds = / fifEdx. (3.2) 

JE JdE Je 

Repeating this for all E G Ph, extending each ip E to zero outside of E, and 
summing in E, yields 

5^1/ { K ( x )^ u - Vip + U(p)dx- / ip(K(x)Vu ■ fj,)ds} 

= J2 I f^ dx - ( 3 - 3 ) 



where 

<p = e {v G L 2 (Q)\v\ E G H\E),VE G P h }. 

E&P h 

Since u is the solution of the problem fll.ip - fll.2p . we know that u satisfies 
the Dirichlet boundary condition on dQ. In addition, it is known that u 



belongs to H(A,£l). From Lemma 13.11 we know that u and (Vu ■ n) are 
weakly continuous across the element interfaces to the variational formulation 
in a weak sense, which yields 

X) { / {K(x)Vu-V<p + wp)dx- (p(K(x)Vu- n)ds} 
' (<p)[K(x)Vu- n}ds + [u](K(x)V<p ■ n)ds 

Tint J Tint 

/ u(K(x)V<p ■ n)ds + f a^-[K(x)Vu ■ n}[K(x)Vip ■ n]ds 



E&P h 

+ 



an Jr int P c 



= [ f<pdx, V<peH 2 (P h ). (3.4) 
Jn 

Combining (12. 9 p and (I2.10p gives 

^2 { J (K(x)Vu-Vip + uip)dx 

(<p(K(x)Vu ■ fj.) - u(K(x)Vip ■ (j,))ds 



E&P h 



dE 

+ f ({cp)[K(x)Vu-n) - (u)[K{x)Vip ■ n))ds 
f h x 

+ a-^[K(x)Vu ■ n][K(x)V(p ■ n]ds, 

/ ftpdx, G H 2 (P h ). 
Jn 



(3.5) 



Applying the density of H 2 (Ph) in V, we complete the proof. □ 
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3.1. Continuity property 

Theorem 3.2. Let B(-, •) be the bilinear form as defined in \2. 7| j. If a > 0, 

then there exists M > such that 

\B(u,v)\ < M\\\u\\\ \\\v\\\, Wu,v e V, (3.6) 
Proof. By the definition of the average, we can obtain 

/ (v)[K(x)Vu ■ n]ds = ~ ^2 [ v[K(x)Vu ■ n]ds, u,v E V. (3.7) 

And by the definition of B(-,-), then we get 
B(u,v) = E{/ ( K ( x )^ u ■ + uv)dx 

E€P h JE 

— I v(K(x)Vu ■ n)ds + / u(K(x)Vv ■ n)ds 

JdE JdE 

+ - I v[K(x)Vu ■ n]ds — — f u[K(x)Vv ■ n]ds\ 
2 JdEnr int ' 2 JoEnV int ' J 

h x f 

+ a— / [K(x)Vu ■ n][K(x)Vv ■ n}ds. (3.8) 
Applying the Schwarz inequality, f)2.3p and f!2.4p . we get 

max |3,— + 



■{ E (ll u ll* + ^il^) Vu -Mll^v 2( a £ )) 

EGP h y 

+ a-\\[K(x)Vu.n]\\l 2{rint) } 

{ E (ll«ll2 + z?ll^(*)v«-^l&-i/. (0 B)) 

h x -\ 1/2 

+a-||[^(x)V,-n]||i 2(rint ; 
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Denote M = max <^3, ^L, + 1 j, where C = max{-^, 1}, hence the proof 
is completed. □ 

3.2. The Inf-sup condition 
3.2.1. The auxiliary problems 

Given an arbitrary u G V, find for every E £ P h the function z E , such 
that 

-V • {K{x)Vz E ) + z e = 0, in E, 

(3.9) 

K(x)Vz E ■ M = K(x)Vu ■ n, on dE. 
The equivalent variational formulation of (13. 9p is that given u G V, find 
z E G H l (E) such that 

(z E ,v)* >E = [ (K(x)Vu ■ rfiovds, veH\E), (3.10) 

JdE 

where (•, •) 3f ^ E denotes the inner product in H 1 ^), and it can be proved that 
(v, f)*,E = \\v \\l easily. By the generalized Lax-Milgram theorem it follows 
that the problem (13.101) has a unique solution z E G H l (E). 

Remark 3.1. Substituting z E and u for v in A3.10\) , we obtain the following 
two identities: 

\\z E \\l = / {K(x)Vu ■ nhoZEds, 

JdE 

{z E: u)* tE = / (K(x)Vu ■ n)>y uds. 

JdE 

From [16|, |20j, we know that the following result holds. 

Theorem 3.3. Given u G V, let z E = z E (u) be the unique solution to Ii3.10\) . 
then the following relation holds: 

\\z E \\* = \\K(x)Vu ■ n\\ H -i/2 {dE) . 
10 



3.2.2. Inf-sup condition on the space V 

In this section, we prove that bilinear form B(-,-) satifies the Inf-Sup 
condition with respect to the norm ||| • |||, defined by f |2.2j) . Let us introduce 
the extension operator ^> e '■ H X (E) — > V, 

x / \ I v E , in E, 

[ 0, in n\E. 

Hence, given a function u G V, we can solve A3. 10[) for a set of functions 
ze(u) and construct a function u G V, such that 

u = u + (3 *e(*e), (3.H) 

where G R. 

Lemma 3.2. Given u G V , then for every G R i/iere exists a strictly 
positive £i = such that 

\\\u\\\ < 

Proof. Substitution of the definition of u into (12. 2p and recalling from (]3.9[) 
that 1T(:e)Vze • n = K(x)Vu ■ on dE, we obtain 



\u\ 



|2 



^ {||m||^ + 2(m,/3^) #)B + \\pz E \ 



+ (1 + /3) 2 ^||K(x)Vu • n\\%- l/2(9E) } 
+ a^(l + /3) 2 ||[K(x)V M -n]||i 2(rmt) . 



Using the Schwarz inequality, triangle inequality, Theorem 13. 3[ and Young's 
inequality as follows: 

2(u,z e )*,e < s\\u\\l + ^\\z E \\l, e > 0, (3.12) 
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we obtain (here taking e = 1) 



\u\ 



< £ {2|kll^+((l + /3) 2 + ^)^||^)V.^||^ 1/2(a , ) } 



|2 ^ 



+ (l + /3)V^||[ir(x)Vn-n]||i 2(rmt) . 



Thus, the assertion holds with 



fc = y / max{2,(l + /3)2 + ^!}. 

□ 

Lemma 3.3. Given u £ V, then there exists £2 = £2(0") ^>p) > sttc/i t/iat 

s(u,u) > 6IIMII 2 - (3.13) 

Proof. By replacing f by m in the definition of B(u,v), and recalling that 
K(x)Vz E ■ M = K(x)Vu ■ /x, we get 

B(u,u) — £ {ikll* + /^( M 5 -z^)*^ - /? / z E (K(x)Vu ■ fx)ds 

+ P u(K(x)Vu- n)ds\ + /3 (z E )[K(x)Vu ■ n]ds 
JdE } Jr int 

-(3 [ (u)[K(x)Vu-n]ds + a^{l + /3)\\[K{x)Vu-n]\\l 2ir . t) . 

•/I.n, P 

For simplicity, the traces JqZe and 70 u have been denoted as and w, 



respectively. Now, using the identities given in Remark 13.11 we can rewrite 
the above expression as follows 



B(u,u) = {\Ml + 2P(u,z E )^ E -P\\z E \\l} 

EeP h 

+ f (z E )[K(x)Vu-n)ds- I (u)[K(x)Vu- n)ds 

J Tint J Tint 



+ + P)\\[K(x)Vu ■ n]\\l Hr . mt) . (3.14) 



12 



If we take a closer look at the terms involving integrals over r int , we see 
that 



I (z E )[K(x)Vu.n]ds>- 1 -^ \M\l-^\\mx)Vu'n]\\h^) 

J^int E£P h 

-0 I (u)[K(x)Vu-n}ds>- 1 -^ Yl ll^ll*- Y ll[ ^ (x)VM ' n]ll W int )- 

•'Tint T?r- EX 



E&P h 

Now, back substitution of these two results into (13.141) . yields 

B(u,u) > ^{(i-M)||„||2 + 2/3(tt,^-(/3+M)|M2} 

+ (a^(l + /3)-|/3|)||[ir(x)Vn-n]||i 2(rint) . 
Using (I3.12p . Cauchy-Schwarz inequality and Theorem 13.31 we obtain 
B{u,u) > J2 \{l-e\/3\-^)\\ur 
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- (P+ l ^+ l ^)\\K(x)Vu. »\\l_ 1/2{9El 
+ (^(l + /3)-|/3|)||[^)V«-n]||| 2(rint) . 

With (5 < 0, we have 



B(u,u) > 



EeP h 

+ (^ 1 -^)ll^)V«-Mll^-v^)} 
+ (^(l-|/3|)-|/3|)ll[^)V W -n]||i 2(rmt) . 

The second term in the right hand side is only positive for e > 4/3. If we 
take e = 2, then we get 



4 



\u\ 
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It is clear that, given the parameters a, A, £, v and 9, we can always find a 
coefficient (3 such that there exists a £2(0", h,p), denoted by 

£ 2 = mm |1 - — , 1 - - j, (3.15) 

that satisfies the inequality A3. 13|) . □ 

Theorem 3.4. Given a > 0, there exists 7 = 7(0", h,p) > such that 

sup ^n^TTr ^ 7lNII> VuGV. (3.16) 
uev\{o} 1 1 Pill 

Proof. By definition of the supremum, we can obtain 

sup M>^M, v» e V. (3.17) 

nev\{o} 1 1 p||| || pill 

where u is defined by (13. lip . Next, by applying Lemmas 13.21 and I3.3[ we 
obtain 

-P ^^>^^>F^4lll«lll. VueV. (3.18) 
»ev\{°} 1 1 Pill 1 1 Pill £i{(T,h,p) 

Taking 7 = £ 2 /£i, we finish the proof. □ 

Corollary 3.1. If\ = u = 9 = ( = 0, then the inf-sup coefficient 7 is a 
constant. 

Proof. For simplicity, we set a = 1. Choosing = 4/10 and using (j3.15|) . 



we have £1 = V288/10 and £ 2 = 1/10, respectively, and it follows that 



7 = 1/V288. □ 

Corollary 3.2. If A = v and 9 = (, then for h x /p^ < 1 the coefficient 7 is 
bounded a constant C > 0. 



14 



Proof. We still set that a = 1, but now we choose (3 = 4h u /10p 8 . If we take 
h x jpf < 1, we obtain the following inequalities from (13 . 1 5[) 



V228 1 

6 -^o"' 6 -To" 



Hence, we conclude that 7 > 1/V228. □ 
3. 3. Existence and uniqueness 

Lemma 3.4. If f G L 2 (f2) ; i/ien i/iere exists a unique solution w £ V to the 
VBVP / TO)) i/iat zs a solution to the problem [Tl\) - [r^) . 

Proof. First, we introduce the classical variational formulation of the model 
problem fO]) - (Q]) in 

Find w G #o(ft), s.t., A(w, v) = L{v), Vf G #o(0), (3.19) 

where L(v) is defined by (12. 8j) and the bilinear form y4(u>,t>) : if 1 (fi) x 
H l (Q) ->■ R is defined by 



i>) 



/ Vw • Vf + wv)dx. 

Jn 



By the generalized Lax-Milgram theorem and by equivalence of this formu- 
lation to the problem (ll.ip -( |L2|) . we know that if / G L 2 (Q) there exists a 
unique solution w G Hq(Q) fl H(A, Q) C V to (13.191) that satisfies the model 



problem in a distributional sense. Consequently, by Theorem I3.1[ we know 
that w G V is a solution to the VBVP (I2.6P as well. 

Thus, following the existence theory for the continuous variational formu- 
lation, we can easily prove the existence of a solution to the VBVP (12. 6p , and 
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we omit the detail of the proof. Also, we know that the solution is unique 
because the bilinear form B(u,v) of (12.71) is positive definite, i.e., 



B(v,v)= N* + ^ll[^) V ^^]llW int )> ' V»6V\{0}. 



3.4- Stability 

The last requirement to ensure well posedness of the weak formulation 
is stability, i.e., the VBVP (12. 6 p is continuously dependent of the solution 
on the input data. Corresponding to Corollary 13.11 and 13.21 we have two 
propositions 13.11 and 13.21 as follows. 

Proposition 3.1. If a = 1, A = ( = 0, and v = 9 = 0, then the solution to 
the VBVP 112.6}) depends continuously on the input data, i.e., given a small 
perturbation 5f G L 2 (Q), then there exists a unique perturbation 5u G V such 
that 



where C is a constant independent of h and p. 

Proof. Let 5f G L 2 (Q) be a perturbation in the input data /. Consequently, 
since the problem under consideration is linear, this leads to a perturbation 
5u G V in the solution u, which satisfies 



E&P h 



□ 



||HI^)<IIIHII<^II^IU 2( n), 




(3.20) 



Using Theorem I3.4[ we have 




— sup 

7 dgv\{o} 




(3.21) 
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Using f)3.20p . (13.211) and the Cauchy-Schwarz inequality, we obtain 

IIHII<- sup <1||^|| (322) 

7 »ev\{0} III^IH 7 

lia — 1,X — u — 9 — ^ — 0, using Lemmas 13.21 and 13. 3[ Theorem I3.4[ 



Corollary 13. II and f)3.22p . we can take C = - = \/228 to complete the proof. 

□ 

Proposition 3.2. If A = v, 9 = (, and h x /p^ < l,then the solution u to the 
VBVP 112.6}) depends continuously on the input data, i.e., given a perturbation 
Sf G L 2 (Q), there exists a unique perturbation Su G V such that 

\\Su\\ HHPh) <\\\5u\\\<C\\5f\\ LHnh 

where C is a constant independent of h and p. 



Proof. Following the proof of Proposition 13.11 and using Corollary 13.21 we 

get 

\\Su\\H H P h) < hsfhHny (3.23) 
Therefore, the proof is completed. □ 

4. Conclusion 

We propose a new variational formulation by a new discontinuous Galerkin 
technique for a two-dimensional reaction-diffusion problem with Dirichlet 
boundary conditions. Our preliminary study discovers that there is a strong 
and intimate connection between the new variational formulation and dis- 
continuous Galerkin methods for the reaction-diffusion problem with discon- 
tinuous coefficient, we intend to further explore this relation, in particular, 
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to make use of this relation to design better hybrid numerical methods which 
hopefully will not only strongly stable in spatial variable and absolutely sta- 
ble in temporal variable but also be optimally convergent. 
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